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■ ABSTRACT 

, We studied the complete randomness of the angular distribution of gamma- 

^] . ray bursts (GRBs) detected by BATSE. Since GRBs seem to be a mixture 

■ of objects of different physical nature we divided the BATSE sample into 5 
00 ' 

o 



subsamples (shortl, short2, intermediate, longl, long2) based on their dura- 



| tions and peak fluxes and studied the angular distributions separately. We 

. used three methods, Voronoi tesselation, minimal spanning tree and multi- 

■ fractal spectra to search for non-randomness in the subsamples. To investigate 

the eventual non-randomness in the subsamples we defined 13 test-variables 
(9 from the Voronoi tesselation, 3 from the minimal spanning tree and one 
from the multifractal spectrum). Assuming that the point patterns obtained 
from the BATSE subsamples are fully random we made Monte Carlo sim- 
ulations taking into account the BATSE 's sky-exposure function. The MC 
simulations enabled us to test the null hypothesis i.e. that the angular dis- 
tributions are fully random. We tested the randomness by binomial test and 
introducing squared Euclidean distances in the parameter space of the test- 
variables. We concluded that the shortl, short2 groups deviate significantly 
(99.90%, 99.98%) from the fully randomness in the distribution of the squared 
Euclidean distances but it is not the case at the long samples. At the inter- 
mediate group the squared Euclidean distances also give significant deviation 
(98.51%). 
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1 INTRODUCTION 



Recentl y, there is no doubt about the cosmo 



GRBs) flZhang k Meszaros 



2004 



Fox et al. 



ogical origin of the gam ma-ray bursts (hereafter 



200^ 



Meszaros! 



20061 ). Then, assuming a large 



scale isotropy for the Universe, one expects the same property for the GRBs as well. Another 
property, which is also expected to occur that GRBs should appear fully randomly, i.e. 
if a burst is observed it does not give any information about the place of the next one. 
If both properties are fulfilled, then the distribu tion is called complete ly random (for the 
astronomical context of spatial point processes see lPasztor fc Tothl (119951 ) ) . There are several 
tests for checking the complete randomness of point patterns, however, these procedures do 
not always give information for both properties simultaneously. 



There are increasing evidence that all the GRBs do not represent a phys i cally homoge- 
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group ( 


"Couveliotou et al. 


1993; 


Horvath 


2002 


Balazs et al. 


2003 



Horvath 


1998 


Mu 


dieriee et al. 


1998 




iakkila et al. 


Hakkila et al. 


2003 


Horvath et al. 


2ooe 


i). Hence, it is 



worth investigating that the physically different sub groups are also dif f erent in their angular 



distrib utions. In the last years the authors provided (IBalazs et al. 



2000a 



1998 



1999 



Meszaros et al. 



pj) several different tests probing t 



le intrinsic isotropy in the angular sky-distribution 



of GRBs collected in BATSE Catalog (IMeegan et al 



20001 ). Shortly summarizing the re- 



sults of these studies one may conclude: A. The long subgroup (T 90 > 10 s) seems to be 
distributed isotropically; B. The intermediate subgroup (2 s ^ T 90 ^ 10 s) is distributed 
anisotropically on the ~ (96 — 97)% significance level; C. For the short subgroup (2 s > T 90 ) 
the assumption of isotropy is rejected only on the 92% significance level; D. The long and 
the short subclasses, respectively, are distributed diffe rently on the 99.3% significance level. 
(About the definition of subclasses see iHorvathJ (119981 ): Tan is the duration o f a GRB, during 
which time the 90% of the radiated energy is received ( IMeegan et al.ll2000l ).) 



Independently and by different tests, 



Litvin et al. 



(120011 ) confirmed the results A., B. and 



C. with one essential difference: for the intermediate subclass a much higher - namely 99.oy/o 
- significance level of anisotropy is claimed. Again, the short subgroup is found to be "sus- 
picious", but only the ~ (85 — 95)% significance level is r eache d. The long subclass seems to 



be distributed isotropically (but see Meszaros fc Stocekl (120031 )). 



Magliocchetti et al. 



(120031) 



found significant an gular correlation on the 2° — 5° scale for GRBs with Tg < 2s durations. 



Tanvir et al. 



(120051 ) reported a correlation between the locations of previously observed short 
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bursts and the positions of galaxies in the local Universe, indicating that between 10 and 25 
per cent of short GRBs originate at low redshifts (z < 0.025). 

It is a reasonable requirement to continue these tests using more sophisticated procedures 
in order to see whether the angular distribution of GRBs is completely random or has some 
sort of regularity. This is the subject of this article. New tests will be presented here. Mainly 
the clarification of the short subgroup's behaviour is expected from these tests. In this 
paper, similarly to the previous studies, the intrinsic randomness is tested; this means that 
the non-uniform sky-exposure function of BATSE instrument is eliminated. 

The paper is organized as follows. In Section [2] the three new tests are described. This 
Section does not contain new results, but - because the methods are not widely familiar 
- this minimal survey may be useful. Section [3] contains the statistical tests on the data. 
Section @] summarizes the results of the statistical tests, and Section [5] presents the main 
conclusions of the paper. 



2 MATHEMATICAL SUMMARY 
2.1 Voronoi tesselation (VT) 

The Voronoi diagram - also known as Dirichlet tesselation or Thiessen polygons - is a funda- 



ment a 



structure in computational geometry and arises naturally in many different applica- 



tions (1Voronoilll908l : IStoyan fc Stoyanlll994l ). Generally, this diagram provides a partition of 
a point pattern ("point field", also "point process") according to its spatial structure, which 
can be used for analyzing the underlying point process. 

Assume that there are N points (N ^> 1) scattered on a sphere sur face with an unit 



radiu s. One says that a point field is given on the sphere. The Voronoi cell (IStoyan fc Stoyan 



1994] ) of a point is the region of the sphere surface consisting of points which are closer to 
this given point than to any other ones of the sphere. This cell forms a polygon on this 
sphere. Every such cell has its area (A) given in steradian, perimeter (P) given by the length 
of boundary (one great circle of the boundary curve is called also as "chord"), number of 
vertices (N v ) given by an integer positive number, and by the inner angles (c^; i — 1, N v ). 
This method is completely non-parametric, and therefore may be sensitive for various point 
pattern structures in the different subclasses of GRBs. 

Note that the behaviour of this tesselation method on the sphere surface is quite different 
from that on the infinite plane. This follows from the fact that the polygon's area will not be 
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Short! GRBs (T 90 < 2s, 0.65 < P 25 6 < 2) 




Figure 1. Voronoi tesselation of the short GRBs (Shortl sample) in the 0.65 < P256 < 2.00 peak-flux range in Galactic 
coordinates. The peak-flux is given in dimesion photon/ (cm 2 s). 



independent from each other, because the total surface of the sphere is fixed in An steradian. 
Hence, the spherical Voronoi tesselation is not effected by border effects, and the Voronoi 
diagram becomes a closed set of convex polygons. 

The points on sphere may be distributed completely randomly or non-randomly; the 
non-random distribution may have different cha racters (clust ering, filaments, etc.; for the 



Diggld dl983h). 



survey of these non-random behaviours see, e.g. 

Random and some regular patterns have distributions of one characteristic maxima (uni- 
modal) but with different variances. The multimodality means different characteristic max- 
ima indicating hierarchical (cluster) structure, the number of modes is determined by the 



and to describe its form (for more details see 


Stovan & Stovan (1994 


) and for the as- 


tronomical context 


Coles & Barrow 


(1990): ColeJ (1991 


): Icke & van de Weveaer 


t (1991): 


Ikeuc 


li & Turner ( 


1991): S 


>ubba Rao & Szalav ( 


1992^; 


van de Weveaer 


(1994): 1 


2aninetti 


(1995 


): Doroshkevich et al. 


(1997): 


Yahaei et al. 


(1999); 


Ramella et al. ( 


2001)). 



2.2 Minimal spanning tree (MST) 

Contrary to VT, this method considers the distances (edges) among the points (vertices). 
Clearly, there are N(N — 1)/2 distances among N points. A spanning tree is a system of lines 
connecting all the points without any loops. The minimal spanning tree (MST) is a system of 
connecting lines, where the sum of the lengths is minimal among all the possible connections 
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Shortl GRBs (T 90 < 2s, 0,65 < P 25 6 < 2) 




Figure 2. The MST for the sample in FigQ] 



Kruskal 


1956; 


Prim 


1957) 



19571 ). In this paper the spherical version of MSF 



is used following the original Prim's paper. 

The N — 1 separate connecting lines (edges) together define the minimal spanning tree. 
The statistics of the lengths and the olmst angles between the edges at the vertices can 



cosmology for studying the statistical properties of galaxy samples ( 
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Doroshkevich fc Turchaninov 



200 lh . 



2.3 Multifractal spectrum 

Let denote P{e) the probability for finding a point in an area of e radius. If P(e) scales as 
e a (i.e. P(e) oc e a ), then a is called the local fractal dimension (e.g. a = 2 for a completely 
random process on the plane). In the case of a monofractal a is independent on the position. 
A multifractal (MFR) on a point process can be defin ed as unification of the subsets of 
different (fractal) dimensions (jPaladin fc Vurpianilll9871 ). One usually denotes with f(a) the 
fractal dimension of the subset of points at which the local fractal dimension is in the interval 
of a, a+da. The contribution of these subsets to the whole pattern is not necessarily equally 
weighted, practically it depends on the relative abundances of subsets. The f(a) functional 
relationship between the fractal dimension of subsets and the corresponding local fractal 
dimension is called the MFR or Hausdorff spectrum. 

In the vicinity of i-th point (i = 1,2, ...,N) one can measure from the neighbourhood 
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Multifroctal spectra 
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Figure 3. MFR spectra of simulated (dot-dashed), Longl (dashed), Short 1 (dotted) and Short2 (three-dot-dashed) samples. 
Boxes represent the error of spectrum points derived from Monte Carlo simulations. Note the shift of the maximum of the 
spectrum of the Shortl sample towards higher values in comparison to a = 2, corresponding to the completely random 2D 
Euclidean case. 

structure a local dimension ai ("Renyi dimension"). This measure approximates the dimen- 
sion of the embedding subset, giving a possibility to c onstruct the MFR spectru m which 



characterizes the whole pattern (for more details see iPaladin &: Vulpianil (119870 ). If the 
maximum of this convex spectrum is equal to the Euclidean dimension of the space, then 
in classical sense the pattern is not a fractal, but the spectrum remains sensitive to the 
non-randomness of the point set. 
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3 STATISTICAL TESTS ON THE DATA 

The three procedures outlined in Section [2] enable us to derive several stochastic quantities 
well suited for testing the non-randomness of the underlying point patterns. 

3.1 Input data and the definition of samples 

Up to the present the most comprehensive all-sky survey of GRBs was done by the BATSE 
experiment on board of the CGRO satellite in the period of 1991-2000. In this period the 
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Table 1. Tested samples of BATSE GRBs. 



Sample Duration Peak flux Number 

[s] [photons cm~ 2 s —1 ] of GRBs 



Short 1 T 90 < 2 s 

Short2 T 90 < 2 s 

Intermediate 2 s < Tgo < 10 s 

Longl T 90 > 2 s 

Long2 T 90 > 10 s 



0.65 < P 25 6 < 2 261 

0.65 < P 256 406 

0.65 < P 256 253 

0.65 < P 256 < 2 676 

0.65 < P 256 966 



experiment collec ted 2704 well justifie d burst events and the data are available in the Current 
BATSE Catalog (IMeegan et al.lboooh 



Since there are increasing evidence (IHorvath et al.l (120061 ) and references therein) that the 



GRB population is actually a mixture of astrophysically different phenomena, we divided 
the GRBs into three groups: short (T go < 2s), intermediate (2s ^ T go ^ 10s) and long 
(T 90 > 10s). To avoid the problems with the changing detection threshold we omitted the 
GRBs having a p eak fl ux P 2 56 ^ 0.65 photons cm~ 2 s _1 . This truncation was proposed by 



Pendleton et al. 



(119971 ). The bursts may emerge at very different distances in the line of 



sight and it may happen that the stochastic structure of the angular distribution depends 
on it. Therefore, we also made tests on the bursts with P256 < 2 photons cm~ 2 s -1 in the 
short and long population, separately Table [1] defines the 5 samples to be studied here. 



3.2 Definition of the test-variables 



The randomness of the point field on the sphere can be tested with respect to different 
criteria. Since different non-random behaviours are sensitive for different types of criteria 
of non-randomness, it is not necessary that all possible tests using different measures reject 
the assumption of randomness. In the following we defined several test-variables which are 
sensitive to different sto chastic properties of the underlying point pattern, as proposed by 
Jl998f ). 



Wallet fc Dussert 



3.2.1 Voronoi tesselation 

Any of the four quantities characterizing the Voronoi cell, i.e. the area, the perimeter, the 
number of vertices, and the inner angles can be used as test- variables or even some of their 
combinations, too. We defined the following test-variables: 

- Cell area A; 

- Cell vertex (edge) N v ; 
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- Cell chords C; 

- Inner angle af, 

- Round factor (RF) average RF av = 4nA/P; 

- Round factor (RF) homogeneity 1 — a ^ av ^ 

- Shape factor A/P 2 ; 

- Modal factor a(aii)/N v \ 

- The so-called "AD factor" defined as AD = 1 — (1 — a(A)/(A})-\ where a (A) is the 
dispersion and (A) is the average of A. 



3.2.2 Minimal spanning tree 

To characterize the stochastic properties of a point patters we use three quantities obtained 
from a MST: 

- Variance of the MST edge-length (t{Lmst)] 

- Mean MST edge-length Lmst\ 

- Mean angle between edges cvmst- 



3.2.3 Multifractal spectrum 

Here the only used variable is the f{a) multifractal spectrum, which is is a sensitive tool for 
testing the non-randomness of a point pattern. 

An important problem is to study the sensitivity (discriminant power) of the different 
parameters to the different kind of regularity inherent in the point pattern. In the case of 
a fully regular mesh, e.g., A is constant and so AD = 0, a{a,i) = and both are increasing 
towards a fully random distribution. In case of a patchy pattern the distribution of the 
area of the Voronoi cells and the edge distribution of MST become bimodal reflecting the 
average area and edge length within and between the clusters, in comparison to the fully 
random case. In a filamentary distribution the shape of the areas becomes strongly distorted 
refl ecting in an increase of t he modal factor in comparison to the case of patches. 



Wallet fc Dussert 



(119971 ) investigated the power of the Voronoi tesselation and minimal 
spanning tree in discriminating between distributions having big and small clusters, full ran- 
domness and hard cores (random distributions but the mutual distances of the points are 
constrained by the size of a hard core), respectively. They concluded that Voronoi round- 
ness factor did not separate small clusters and hardcore distributions, and roundness factor 
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homogeneity did not distinguish between small clusters and random distributions, nor be- 
tween random and hardcore distributions. MST has a very good discriminant power even in 
the case of hardcore distributions with close minimal interpoint distances. 

Since the sensitivity of the variables are different on changing the regularity properties of 
the underlying point patterns one may measure significant differences in one parameter but 
not in the other even in the case when these are correlated otherwise. It is not a trivial issue. 
In most cases one needs extended numerical simulations to study the statistical significance 
of the different parameters. 



3.3 Estimation of the significance 

Let us denote with £ one of the thirteen test- variables defined in Section 1X21 The probability 
that £ < x occurs is given by P(£ < x) = F(x), where F(x) is the probability distribution 
function. We approximated F(x) numerically by the F n (x) empirical probability distribution 
function which can be calculated by F n (x) = k/n where n is the number of simulations and 
k is the number of cases for which the simulated £ < x holds. 

Similarly, the probability that £ is within the interval [x\ \ X2] can be obtained by F n (xi) — 
F n (x2) = (x 2 — Xi)/n. Then the (3 probability that £ is outside this region is given by 
(3 = 1 — (x2 — Xi)/n. In the following we suppose that the [aJi;^] interval is symmetric to 
the x sample mean. To obtain the empirical distributions of the test-variables we made 200 
simulations for each of the five samples. The number of the simulated points were identical 
with those of the samples defined in Section 13.11 

We generated the fully random catalogs by Monte Carlo (MC) simulations of fully ran- 
dom GRB celest i al po s itions and tak i ng int o account the BATSE sky-exposure function 



flFishman et al. I (119941 ) . iMeegan et all ((20001)). 

Assuming that the point patterns obtained from the five samples, defined in Table 1, are 
fully random we calculated the probabilities for all the 13 test-variables selected in Section 
13.21 Based on the simulated distributions we computed the level of significance for all the 
13 test- variables and in all samples defined. 
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4 DISCUSSION OF THE STATISTICAL PROPERTIES 
4.1 Significance of independent multiple tests 

In Section 13.31 we calculated numerical significance for the tests assuming they were per- 
formed individually. The calculated significance levels are given in Table O In the reality, 
however, these figures would mean significance at a certain level if one performed only that 
single test. Assuming that all the single tests were independent the P n (m) probability that 
among n trials at least m will resulted significance only by chance at a certain level is given 
by the following equality: 

n 

Pn(m)=Y,P£, (!) 

k=m 

where PJ} is the binomial distribution giving the probability of k successes among n trials: 

In the equation given above p means the probability that a single test has given significant 
result only by chance. It is easy to see that this equation resulted in -P n (l) = 1 — (1— p) n ~ np 
which gives a significance of 1 — np approximately instead of 1 —p. It means e.g. that a single 
test resulted in 1 — p = 0.95 significance is reduced to 1 — 0.95 13 = 0.49 if one performed 
n = 13 independent tests but only one resulted in 1 — p = 0.95 significance. 

Inspecting Table [2] listing the calculated numerical significance of single tests one can 
infer that shortl sample has 4 tests with significance of 1 — p = 0.95 ore more. Taking 
into account the calculations at the end of the previous paragraph, however, we have to 
emphasize that the calculated numerical significance, based on the individual probability 
distribution of the test- variables separately, does not have its original meaning. Significance 
refers to the certainty rejecting the null hypothesis on the basis of the bunch of tests as a 
whole. Applying Equation ([1]) with m = 4 and n = 13 one gets a significance of 99.69 %. 
Applying the same sequence of arguments to the shortl sample one may get a figure of only 
86.46 % (m = 2 and n = 13). In the case of the intermediate, longl and long2 samples one 
can not get figures above the 95 % significance level. 

One may have a serious concern, however, with the results obtained above. Namely, the 
basic requirement of the independence of the single tests is not fulfilled in our case. In the 
contrary, there are strong correlations between the test-variables in Table El In the next 
subsection we try to outline an approach which takes into account the correlations between 
the test-variables and overcomes this difficulty 
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Table 2. Calculated significance levels for the 13 test-variables and the five samples. A calculated numerical significance 
greater than 95% is put in bold face. 

Name var shortl short2 interm. longl long2 

Cell area A 36.82 29.85 94.53 79.60 82.59 

Cell vertex (edge) N v 36.82 87.06 2.99 26.87 7.96 

CeU chords C 47.26 52.24 18.91 84.58 54.23 

Inner angle at 96.52 21.39 87.56 37.81 63.18 

RF average A-kA/P 65.17 99.98 33.83 10.95 86.07 

RF homogeneity 1 - tT( J f/°" ) 19.90 24.38 58.71 55.72 32.84 

Shape factor A/P 2 °" 91.04 94.03 90.05 55.22 63.68 

Modal factor a(ai)/N v 97.51 1.99 7.46 56.22 8.96 

AD factor l-(l-^^) _1 32.84 25.37 11.44 95.52 52.74 

MST variance v(L MST ) 52.74 38.31 22.39 13.93 59.70 

MST average L M st 97.51 7.46 89.05 56.72 8.96 

MST angle a MST 85.07 14.43 36.82 73.63 60.70 

MFR spectra f(a) 95.52 96.02 98.01 73.63 36.32 



Binomial test 


(Eq. QJ with p = 0.05) 


99.69 


86.46 


51.33 


51.33 


Squared Euclidean 


distance 


99.90 


99.98 


98.51 


93.03 36.81 



4.2 Evaluation of the joint significance levels 

We assigned to every MC simulated sample 13 values of the test variables and, consequently, 
a point in the 13D parameter space. Completing 200 simulations in all of the subsamples 
we get in this way a 13D sample representing the joint probability distribution of the 13 
test-variables. Using a suitable chosen measure of distance of the points from the sample 
mean we can get a stochastic variable characterizing the deviation of the simulated points 
from the mean only by chance. An obvious choice would be the squared Euclidean distance. 

In case of a Gaussian distribution with unit variances and without correlations this 
would resulted in a \ 2 distribution of 13 degree of freedom. The test-variables in our case 
are correlated and have different scales. Before computing squared Euclidean distances we 
transformed the test- variables into non-correlated ones with unit variances. Due to the strong 
correlation between some of the test-variables we may assume that the observed quantities 
can be represented with non-correlated variables of less in number. Factor analysis (FA) is 
a suitable way to represent the correlated observed variables with non-correlated variables 
of less in number. 



Since our test- variables are stochastically dependent following IWallet fc Dussertl (119981 ) 



we attempted to represent them by fewer non-correlated hidden variables supposing that 

k 

Xi = Y^ <kjfj + ^ i = 1, • • • , V ; k < p . (3) 

j'=i 

In the above equation Xi,fj,Si mean the test- variables (p = 13 in our case), the hidden 
variables and a noise-term, respectively. Equation ([3]) represents the basic model of FA. 
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After making some reasonable assumptions (IKendall fc Stuartlll973l ). k can be constrained 
by the following inequality: 



fc^(2p + l-y8p + l)/2 (4) 

which gives k ^ 8.377 in our case. 

Factor analysis is a common ingredient of professional statistical software packages 
(BMDP, SAS, S-plus, SPSsQ, etc). The default solution for identifying the factor model is to 
perform principal component analysis (PCA). We kept as many factors as were meaningful 
with respect to Equation (j3J). Taking into account the constraint imposed by Equation 
we retained 8 factors. In this way we projected the joint distribution of the test- variables in 
the 13D parameter space into a 8D one defined by the non-correlated fa hidden variables. 

The f) hidden variables in Equation (j3J) are non-correlated and have zero means and 
unit standard deviations. Using these variables we defined the following squared Euclidean 
distance from the sample mean: 

d 2 = + /! + ... + fl (k = 8 in our case) . (5) 

If the fj variables had strictly Gaussian distributions Equation (jSJ) would define a x 2 variable 
of k degrees of freedom. 



4.3 Statistical results and their interpretations 

In addition to the significance obtained by the binomial test in Subsection 14.11 using the 
distribution of the squared Euclidean distances, defined by Equation (jHJ), one can get further 
information whether a BATSE sample represented by a point in the parameter space of 
the test-variables deviates only by chance or it significantly differs from the fully random 
distribution. 

In all categories (shortl, short2, intermediate, longl, long2) we made 200, altogether 
1000, simulations. We calculated the d 2 squared distances for all simulations and compared 
them with those of the BATSE samples in Table [TJ Figure H] shows a histogram of the 
simulated squared distances along with those of the BATSE samples. Full line represent a 
X 2 distribution of k = 8 degree of freedom. Figure H] clearly shows that the departures of 
samples shortl and short2 exceed all those of the simulated points. The probabilities, that 
these deviations are non-random, equal 99.9% and 99.98%, respectively. 



1 BMDP, SAS, S-plus, SPSS are registered trademarks 
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Figure 4. Distribution of the squared Euclidean distances of the simulated samples from the stochastic mean of the /; hidden 
variables (factors) in the 8D parameter space. There are altogether 1000 simulated points. Full line marks a \ 2 distribution of 
8 degree of freedom, normalized to the sample size. The distances of the BATSE samples are also indicated. The departures of 
samples " shortl" and " short2" exceed all those of the simulated points. The probabilities, that these deviations are non-random, 
equal 99.9% and 99.98%, respectively. 



The full randomness of the angular distribution of the long GRBs, in contrast to the 
regularity of the short and in some extent to the intermediate ones, points towards the 
differences in the angular distribution of their progenitors. The recent discovery of the af- 
terglow in s ome short GRB s indicates that these events are associated with the old stellar 
population (IFox et al.ll2005l ) accounted probably for the mergers of compact binaries, in con- 
trast to the long bursts resulting from the collapse of very massive stellar objects in young 
star forming regions. The differences in progenitors reflects also the differences between the 
energy released by the short and long GRBs. 

Unfortunately, little can be s a id on the physical nature of the intermediate class. The 
statistical studies (jliorvath et al.l ^OOfj ) and the references therein) suggest the existence of 
this subgroup - at least from the purely statistical point of view. Als o the non-random sk y 
distribution is occurring here. But its physical origin is fully open yet (IHorvath et al.ll200q ). 
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5 SUMMARY AND CONCLUSIONS 

We made additional studies on the degree of the randomness in the angular distribution of 
samples selected from the BATSE Catalog. According to the T 90 durations and P 256 peak 
fluxes of the GRBs in the Catalog we defined five groups: shortl (T 90 < 2 s & 0.65 < P 256 < 
2), shortl (T 90 < 2 s & 0.65 < P 256 ), intermediate (2 s ^ T 90 ^ 10 s & 0.65 < P 256 ), longl 
(T 90 > 2 s & 0.65 < P 256 < 2) and longl (T 90 > 10 s k 0.65 < P 256 ). 

To characterize the statistical properties of the point patterns, given by the samples, 
we defined 13 test- variables based on the Voronoi tesselation (VT), Minimal spanning tree 
(MST) and Multifractal spectra. For all five GRB samples defined we made 200 numerical 
simulations assuming fully random angular distribution and taking into account the BATSE 
exposure function. The numerical simulations enabled us to define empirical probabilities for 
testing the null hypothesis, i.e. the assumption that the angular distributions of the BATSE 
samples are fully random. 

Since we performed 13 single tests simultaneously on each subsamples the significance 
obtained by calculating it separately for each test can not be treated as a true indication 
for deviating from the fully random case. At first we supposed that the test-variables were 
independent and making use the binomial distribution computed the probability of obtaining 
significant deviation in at least one of the variables only by chance. In fact, some of the test- 
variables are strongly correlated. To concentrate the information on the non-randomness 
experienced by the test-variables, we assumed that they can be represented as a linear 
combination of non-correlated hidden factors of less in number. Actually, we estimated k = 8 
as the number of hidden factors. Making use the hidden factors we computed the distribution 
of the squared Euclidean distances from the mean of the simulated variables. Comparing the 
distribution of the squared Euclidean distances of the simulated with the BATSE samples we 
concluded that the shortl, short2 groups deviate significantly (99.90%, 99.98%) from the fully 
randomness, but it is not the case at the long samples. At the intermediate group squared 
Euclidean distances also give significant deviation (98.51%). 
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